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DELAY  ANALYSIS  OF  INTERACTING  QUEUES  WITH  AN  APPROXIMATE  MODEL 

*  by 

Anthony  Ephremides  and  Rong-Zhu  Zhu 

The  University  of  Maryland  China  Signals  &  Communications  Co. 

ABSTRACT 

An  approximate  model  of  coupled  Markov  chains  is  proposed  and  analyzed  for  a 
slotted  ALOHA  system  with  a  finite  number  of  buffered  nodes.  This  model  differs  from 
earlier  ones  in  that  it  attempts  to  capture  the  interdependence  between  the  nodes. 

The  analytical  results  lead  to  a  set  of  equations  that,  when  solved  numerically, 
yield  the  average  packet  delay.  Comparison  between  computational  and  simulation 
results  for  a  small  number  of  nodes  show  excellent  agreement  for  most  throughput 
values,  except  for  values  near  saturation.  Numerical  comparisons  for  a  two-node 
system  show  that  a  non-symmetric  loading  of  the  system  provides  better  delay- 
throughput  performance  than  a  symmetric  one. 

1.  INTRODUCTION 

Numerous  papers  have  appeared  in  the  literature  that  study  and  analyze  the  pro¬ 
tocol  of  slotted  ALOHA.  However  most  of  them  either  concentrate  on  capacity  analysis 
or  assume  an  infinite  user  population  that  implies  no  queueing  of  packets  at  the  ter¬ 
minals.  The  case  of  a  finite  number  of  terminals  that  receive  packets  from  exogenous 
sources  and  that  provide  buffer  space  for  queueing  is  of  interest  because  it  is  both 
practical  and  mathematically  leading  to  a  challenging  problem.  In  this  paper  we  con¬ 
sider  this  problem.  ,t 

In  [l]  an  approximate  model  was  introduced  that  proved  analytically  tractable 
but  did  not  lead  to  good  agreement  with  simulation  results.  Nevertheless  it  repre¬ 
sented  an  attempt  to  capture  the  interdependence  between  the  nodes,  contrary  to  the 
often  assuned  practice  of  considering  them  statistically  independent,  and  formed  the 

basis  for  similar  approximate  analyses  in  related  problems  [2,3].  Here  we  propose  a 

1 

*This  work  was  supported  in  part  by  ONR  Grant  N00014-84K-614  and  by  NSF  Grant 
85-00108. 


refinement  of  the  model  presented  in  [l  ]  which  tracks  closer  the  Interaction  between 
the  nodes  and  thus  leads  to  the  expectation  of  better  agreement  with  the  actual 
system.  This  expectation  is  confirmed  by  comparison  of  the  analytical  results  to 
those  of  a  simulation  that  was  carried  out  only  for  a  small  number  of  nodes  (two  and 
three,  respectively). 

The  precise  modeling  of  the  system  we  consider  is  straightforward  in  terms  of  a 
finite-dimensional,  infinite,  discrete  random  walk  with  appropriate  boundaries.  It 
is  in  fact  the  presence  of  the  boundaries  that  causes  complications  and  makes  the 
analysis  difficult  [l ,4  ] .  Even  approximate  models  have  proven  quite  difficult  to 
analyze  [5-7].  The  theoretical  basis  for  establishing  this  as  an  interesting  analy¬ 
tical  problem  can  be  found  in  [8],  where  a  closely  related  problem  of  two  coupled 
processors  was  considered  and  solved  exactly.  In  [9-10  ]  other  approximate  analysis 
techniques  were  proposed  and  met  with  limited  success  in  terms  of  predicting  the 
simulation  results.  In  these  same  papers,  as  well  as  in  [ll],  the  related  problem  of 
system  stability  was  also  considered.  In  [l2 ]  and  in  [l3 ]  an  idea  of  stochastic 
system  domination  has  been  introduced  that  provides  stability  bounds  and  approximate 
derivation  of  the  steady-state  distribution  of  an  ergodic  Markov  chain  respectively. 
There  is  a  very  rich  literature  on  problems  that  stem  from  this  central  problem  of 
analysis  of  a  slotted-ALOHA  system  with  a  finite  number  of  buffered  users.  Thus,  one 
more  effort  in  trying  to  "discover"  a  good  approximation,  as  this  paper  attempts  to 
do,  is  perhaps  justified. 

2 .  The  Problem 

Consider  M  terminals  each  of  which  receives  packets  from  exogenous  sources 
according  to  a  Bernoulli  process  and  lines  them  in  a  queue  in  a  buffer  of  infinite 
capacity.  Hie  arrival  rate  at  the  ith  station  is  A^  and  the  arrival  processes  at 
different  stations  are  statistically  independent.  Time  is  slotted  and  it  takes 
exactly  one  slot  to  transmit  one  packet,  or  for  a  packet  to  arrive.  Thus  A^  is  also 

the  probability  of  an  arrival  at  the  Ith  station  in  any  given  slot. 
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The  ith  terminal  attempts  with  probability  p^  to  transmit  the  head-of-the-line 
packet  in  the  queue  (if  the  latter  Is  non-empty)  through  a  common  channel  to  the 
receiver.  The  usual  ALOHA  assumptions  are  made  concerning  success  or  collision. 
Specifically,  since  we  are  considering  discrete  time,  we  are  assuming  that  a  packet 
that  arrives  at  time  t  joins  the  queue  at  the  corresponding  terminal  and  it  is 
possible,  if  the  queue  is  empty  at  that  instant,  that  it  will  be  transmitted  at  that 
same  instant.  Ternary  information  (success,  idle,  collision)  becomes  available  to 
all  terminals  at  the  next  instant  t+1  and  it  is  possible  to  retransmit  or  delete  the 
packet  at  that  new  time  depending  on  whether  it  was  involved  in  a  collision  or  a  suc¬ 
cessful  transmission. 

To  include  several  variants  of  the  ALOHA  protocol  (such  as  immediate  first 
transmission,  delayed  first  transmission,  etc.  (see  [l4]))  we  distinguish  the  status 
of  a  given  terminal  at  the  beginning  of  a  given  slot  into  three  categories,  idle, 
active  and  blocked.  A  terminal  is  idle  if  there  were  no  packets  in  its  buffer  at 
the  end  of  the  preceeding  slot,  it  is  blocked  if  its  queue  is  not  empty  and  the 
latest  attempted  transmission  was  unsuccessful,  and  it  is  active  if  its  queue  is  not 
empty  but  its  most  recent  attempted  transmission  was  successful.  We  then  let 

r1Ai,  if  i  is  idle 

p^  *  s^,  if  i  is  active 

(/ 

qif  if  i  is  blocked 

For  this  system  which  is  schematically  shown  in  Fig.  1,  we  wish  to  calculate 
either  the  joint  probability  distribution  of  the  queue  sizes  {n^ , i- 1 , . . . ^ }  at  steady 
state  or  the  average  values  of  these  queue  sizes  (and,  therefore,  the  average  packet 
delays  using  Little's  formula). 

The  mathematical  model  for  this  system  consists  of  a  2M  dimensional  Markov  chain 
the  states  of  which  are  represented  by 
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S  “  (S^  , . . .  ,S^;  , • • • jN^) 


where 


S1  “  status  of  the  ith  terminal  (ternary  variable  with  values  «  idle,  active,  and 


blocked) 


and 


*  queue  size  at  the  ith  terminal 

The  transition  probabilities  for  this  chain  can  be  obtained  explicitly  and  so 
can  the  steady  state  equations.  However  their  nunber  is  infinite  and,  owing  to  the 
fact  that  the  transition  probabilities  on  the  boundaries  are  different  from  the  ones 
in  the  interior  of  the  state  space,  traditional  analysis  methods  (for  example,  in 
terms  of  moment  generating  functions)  are  inadequate. 

3.  The  Proposed  Approximate  Model 

The  idea  of  simplifying  the  model  without  losing  track  of  the  Interdependence 
between  the  queue  sizes  of  different  terminals,  that  was  proposed  in  [l],  consisted 
of  considering  the  "status-state"  separately  from  the  "queue-size"  state.  Thus  a 

system  state  was  defined  that  consisted  of  the  triplet  (I,A,B)  where  I  denoted  the 

v 

minber  of  idle  terminals,  A  the  nunber  of  active  terminals,  and  B  *  M  -  (I+A)  the 

nunber  of  blocked  terminals.  Separately,  for  each  terminal  a  node  state  was  defined 

consisting  of  the  pair  (N^,  S^)  where  denoted  the  queue  size  of  the  ith  terminal 

✓ 

and  its  status. 

Transitions  between  the  states  of  the  system  state  variable  are  not  Markovian. 
However,  in  [l ]  they  were  considered  Markovian  with  transition  probabilities  equal  to 

*  we  represent  "idle"  by  0,  "active"  by  1,  and  "blocked"  by  2. 
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the  average  probabilities  of  transition,  where  the  averaging  took  place  over  all 
values  of  the  queue  sizes.  Similarly  the  transitions  between  the  states  of  the  node 
state  variable  are  not  Markovian.  Again,  in  [l  ]  they  were  considered  Markovian  with 
the  transition  probabilities  calculated  as  averages  of  the  true  transition  probabili¬ 
ties  over  all  values  of  the  system  status  states.  Ihe  resulting  balance  equations 
for  both  systems  were  coupled  and  could  be  solved  numerically.  As  mentioned  earlier, 
agreement  with  simulation  results  proved  to  be  unsatisfactory  and  this  was  proved  by 
a  simulation  study  done  at  ETH  [l5 ] . 

In  trying  to  determine  the  reason  for  the  limited  success  of  the  model  proposed 
in  [l],  we  observed  that  the  model  failed  to  track  the  status  of  the  M  terminals  one- 
by-one.  Instead  it  considered  only  the  nimbers  of  active,  idle,  and  blocked  ter¬ 
minals  respectively. 

In  this  paper  we  formulate  a  model  that  again  separates  status-states  from  queue- 
size-states,  in  the  spirit  of  the  model  in  [l],  but  expands  the  state  space  of  the 
system  status  variables  in  order  to  track  every  terminal's  status  separately.  Again, 
however,  the  model  is  approximate  in  that  the  transition  probabilities  between  the 
states  of  each  of  the  two  chains  are  the  averages  of  the  true  transition  probabili¬ 
ties  from  the  original  2M-dimensional  Markov  chains. 

Before  we  proceed  with  the  mathematical  details  we  would  like  to  justify  the 
value  of  this  approximation  a  priori.  It  is  clear  that  in  a  slotted  ALOHA  system  the 
probability  of  success  for  an  attempted  transmission  depends  exclusively  on  the  sta¬ 
tus  of  the  other  terminals.  Ihus  for  the  ith  terminal  the  probability  of  success  is 
given  by  # 

0 

p(i)  -  P.  n  (i-p  )  -  P.  n  (i-r  x  )  n  (i-s  )  n  (i-q  ),  (D 

j*i  3  j*L  J  J  j*L  3  i*i  3 

j  cSl  j  #SA  j  eSB 


where  Sj , 
t ively . 


S.,  Sn  represent  the  sets  of  the  idle,  active,  and  blocked  nodes  respec- 
A  B 

If  the  system  is  "symmetric" ,  that  is  the  quantities  r.,s  ,q  ,1  do  not 

c  J  J  J  J 


depend  on  j,  the  success  probability  becomes  a  function  only  of  the  nimbers  of  ter¬ 
minals  in  the  sets  ST ,  S  ,  and  S  .  This  observation  led  to  the  model  developed  in 
[l].  However,  transitions  between  the  states  of  a  process  that  consists  of  these 
numbers  only,  fails  completely  to  "remember’*  the  status  of  a  given  terminal  in  going 
from  slot  to  slot.  Thus,  although  the  success  probability  might  still  depend  on  the 
nunbers  only,  the  transitions  "average-out"  the  effect  of  the  tendency  of  a  given 
terminal  to  undergo,  or  postpone,  a  change  in  its  own  status.  Of  course,  in  addition, 
if  the  system  is  not  symmetric,  it  is  necessary  to  track  the  status  of  each  terminal 
separately  because  each  terminal  contributes  differently  to  the  success  probability 
depending  on  its  status. 

Thus  we  shall  consider  a  status  variable  S  consisting  of  M  ternary  variables, 

Sj ,S^t ... ,S^,  each  of  which  indicates  the  status  of  the  corresponding  terminal. 

Separately,  we  shall  consider  M  queue  length  variables,  each  representing  the  status 
and  queue  length  of  a  given  terminal  independently  of  the  status  or  queue  lengths  of 
the  other  terminals.  We  shall  consider  these  variables  to  evolve  as  Independent 
Markov  chains  (they  are,  in  fact,  neither  Markov,  nor  independent).  The  transition 
probabilities  of  each  chain,  we  shall  see,  will  depend  on  the  steady-state  probabili¬ 
ties  of  the  other  chains  and  thus  will  take  into  account  a  parametric  (but  not  statistical 
coupling  between  them,  even  though  the  model  is  not  exact. 

3.1  The  System-Status  Markov  Chain 

Let  S  *  (S^,...,Sjj)  be  the  M-dimensional  status  variable,  where 

(o,  if  i  is  idle 


Si  «  1, 


2, 


if  i  is  active 
if  1  is  blocked 


Let  P(S)  be  the  joint,  steady-state  probability  distribution  of  this  random  vec¬ 
tor.  In  Fig.  2  the  exact  states  for  M=2  are  shown.  (Mote  that  the  state  (1,1)  is 

not  achievable  since  it  is  Impossible  for  two  nodes  to  have  had  successful 

6 


transmissions  in  the  same  slot).  The  total  number  of  states  achievable  by  this 
vector  is  easily  calculated  to  be 

I  (%  -  2«- v» 

n-0  n  n=0  n  n=0  n 

where  n  indicates  the  number  of  blocked  nodes,  and  where  the  fact  that  not  more  than 
one  node  can  be  active  is  taken  into  account. 

It  is  assumed  that  the  evolution  of  S  is  Markovian  with  probabilities  that  are 
suitably  estimated.  Their  calculation  is  displayed  in  Appendix  1,  where  we  assume 
=  r^,  =  1.  It  is  important  to  note  that  in  that  calculation  two  quantities  that 
characterize  the  node-chain  (not  yet  described)  are  required.  These  are  calculated 
in  Appendix  2  and  are  the  following: 

P^(l  |l)  =  Pr[queue  size  >  l|  node  i  is  active] 

and 

P^(0j2)  =  Pr [queue  size  =  l|  node  i  is  blocked] 

Thus  the  state  equations  for  this  chain  (and  therefore  the  solution  P(S))  depend 
on  the  above  quantities.  As  we  shall  see,  the  above  quantities,  as  obtained  in 
Appendix  2,  depend  in  turn  on  P( S) .  Thus  all  state  equations  for  both  chains  must  be 
solved  simultaneously. 

3.2  The  Queue-Length  Markov  Chain 

✓ 

Let  be  the  total  nunber  of  packets  (the  queue-length)  at  terminal  i  and  let 
T^  be  the  Indicator  of  blocked  or  unblocked  status,  i.e. 

*  note  that  the  queue  size  is  considered  to  Include  the  packets  waiting  in  the  queue 
to  be  served  plus  the  packet  being  served;  that  is,  it  is  equal  to  the  total  number 
of  packets  residing  at  the  terminal. 
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0,  if  i  is  blocked 

1,  if  i  is  unblocked  (i.e.  active  or  idle) 

The  pair  (T^,N^)  constitutes  the  state  of  the  terminal  and  let  »(T^,N^) -denote 
its  steady  state  probability.  This  state  is  assumed  to  evolve  as  a  Markov  chain.  Of 
special  interest  in  our  calculations,  as  shown  in  Appendix  2,  are  the  quantities 
ir(0,0) ,  *(1,0),  *(1,1)  and  the  conditional  moment  generating  functions  G*(z)  and 
G*(z),  where 

i  Ni 

GZ.  (z)  *  Z  *(T  N  )z  ,  T.  “  0,1 
Ai  N^O 

These  quantities,  as  well  as  the  state  transition  probabilities  that  lead  to  their 
determination,  can  be  calculated  in  terms  of  the  following  "average"  success 
probabilities: 


PB(i)  *  Pr [success  |node  i  is  blocked] 

(2) 

PA(i)  «  Pr [success |noide  i  is  active] 

(3) 

Pj(i)  ■  Pr [success jnode  is  is  idle] 

(4) 

where  the  averaging  is  performed  over  the  status  of  the  other  nodes  (and  for  which 
the  state  probabilities  P(S)  of  the  status  states  are  required,  thereby  leading  to 
the  coupling  of  the  two  sets  of  equations). 

In  fact  it  is  shown  in  the  appendix  that,  for  the  case  of  r^  *  s^  *  1  and  q.^  * 
p^  (case  of  immediate  first  transmission),  we  have 


*(0 ,0) 


J^PjCD 


W1)  +  XiPB(i) 


*(1,0) 


*(1,0)  -  — 


XiPB(i)  -  \?kU) 


XtPB(i)  -  Xi(PI(i)  -  PA(i)) 


(5) 


(6) 
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*(1,1) 


*(0,0) 


(7) 


®q(  1)  "  = 


xix1PI(i) 


(8) 


x^gd)  -  xi(PB(i)  -  pa(D) 


Gl(1)  -  xi+  X1T 


XiPB(1)  "  XiPA(i) 


(9) 


XiPB(1)  “  W1)  "  PA(1)> 


where  X  *»  1  -  X  and  P  *  1  -  P. 

Similar  expressions  can  be  obtained  for  the  much  studied  case  of  delayed  first 
transmission  (i.e.  r^  •»  s^  *  ■  p^). 


4 .  The  Average  Success  Probabilities 

Central  to  the  calculations  in  the  appendices  is  the  knowledge  of  the  average 
success  probabilities  as  defined  above  in  Eqs.  (2-4).  These  are  easily  expressed  in 
terms  of  the  system  state  probabilities  P(Sj,...,S^)  and  by  making  use  of  the  simple 
ALOHA  logic  expressed  in  Eq.  (1),  as  follows: 

1)  for  blocked  node  i 

P  (i)  *  Pr(node  i  success  Inode  i  is  blocked) 
o 

Pr(node  i  success  and  node  i  is  blocked) 

Pr(node  i  is  blocked) 


where 


P(node  i  success  and  node  i  is  blocked) 


_62S— 60S- 

p  z  n  (p/S3X  UiJ)P(S 
1Si“2  jtt  3  J 


S<*1 


(1«> 


.. . ,SM) 
0* 


with  the  summation  ranging  over  all  state  values  for  which  node  i  is  blocked  and  no 
other  node  is  active.  The  exponents  6  are  Kronecker  coefficients.  Furthermore, 
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P(node  i  blocked)  -  Z  P( .S^, . . .S^) 

Si~2 

2)  for  active  node  1 

P^(i)  »  P(node  i  success |node  i  Is  active) 

P(node  i  success  and  node  i  is  active) 
P(node  i  is  active) 

where 


P(node  i  success  and  node  i  is  active) 


Z 

Si 

*1 


n  (p  2sJx  osJ)p(s 

i  j«J  J  1 


•»SM> 


(j*i) 


and 


P(node  i  is  active)  *  Z  P(  Sj  , . .  .S^, . .  .Sj^) 

si“l 

3)  for  idle  node  i 


Pj(i)  =  P(node  i  success jnode  i  is  idle) 

P(node  i  success  and  node  i  is  idle) 
P(node  i  is  idle) 


where 


P(node  i  success  and  node  i  is  idle) 


Z 

Si* 


0 


n  (p^Jx^^-pcSj 


Sj 

<3«> 


and 


P(node  i  is  idle)  -=  Z  P(  S.  , . .  .S, , . .  .Su) 

sA-o  1  1  M 
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5.  The  De lay 


The  total  average  delay  incurred  by  a  packet  from  the  moment  of  its  arrival  to 
the  terminal  until  the  time  of  its  successful  transmission  consists  of  three  com¬ 
ponents,  waiting  time  in  queue,  Wq(i),  "service"  time  from  the  instant  it  arrives  at 
the  head-of-the-line  position  in  the  buffer  until  it  is  successfully  transmitted, 
Wg(i) ,  and  the  actual  transmission  time  (which  we  take  to  be  equal  to  one  slot  for 
"local"  environments). 


G0(1) 

The  component  W  (i)  is  approximated  by —tty  which  is,  equivalently,  given  by 
8 

Pr(i  is  blocked)  .  .  .  ,  .... 

- -7VV - .  The  component  W  (i)  is  obtained  from  Little's  result  by 

PBu)  q 


vu 


where  is  the  average  queue  length  (without  considering  the  blocked  head-of-the- 
line  packet  as  part  of  the  queue).  This  quantity  is  calculated  in  Appendix  2  and  is 
given  by 


Xi\PI(i) 


(  VB(i)  "  XiPA(1))(XiPB(i)  "  W1*  "  PAU))) 
Thus,  for  terminal  i,  the  delay  is  given  by 

Li  G0(1) 

°i“  1+  Xi+  PB(i) 


(10) 


and  the  total  weighted  average  system  delay  D  is  given  by 
M 

I  Vi 

i«l 


M 

1  Xi 
i-1  1 


6.  The  Solution  Algorithm  and  Numerical  Results 

As  explained  in  the  preceding  sections  the  state  equations  that  must  be  solved 
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to  yield  the  system  state  probabilities  P( .S^ , . . . .S^)  and  the  state  equations  for 
the  node  state  probabilities  ir(T^,N^)  are  coupled  and  must  be  solved  simultaneously. 

The  coupling,  as  shown  in  the  appendixes,  occurs  in  terms  of  the  boundary  condition 
probabilities  P±(  1  1 1)  and  P^O  |2)  as  well  as  the  average  success  probabilities  Pg(i), 
PT(i),  PA(i).  Given  the  values  of  these  auxiliary  quantities  each  system  of 
equations  can  be  solved  separately.  Thus  due  to  the  non-linear  nature  of  the 
coupling,  it  is  natural  to  proceed  iteratively,  as  shown  in  Figure  3.  The  iterative 
method  used  is  identical  to  the  one  used  in  [l].  Arbitrary  values  for  one  set  of 
state  probabilities  are  chosen  and  used  to  solve  for  the  other  set.  Then  the  values 
obtained  for  the  second  set  are  used  to  recalculate  those  of  the  first  set,  and  so 
on.  Conditions  for  convergence  are  identical  as  in  [l].  Of  course  it  must  be  assu¬ 
med  that  the  values  of  the  arrival  rates  X^  are  such  that  the  system  of  the  queues  is 
stable  (ergodic).  As  conditions  for  ergodicity  are  not  known  except  for  the  case  of 
M»2  [9-13],  we  used  rather  small  values  for  the  X^'s  to  avoid  running  into  the 
unstable  region  in  which  the  nunerical  results  would  be  meaningless.  In  fact,  the 
region  near  saturation  is  expected  to  show  poor  agreement  with  simulation  exactly 
because  the  stability  boundary  (which  is  not  precisely  known)  may  be  crossed  in  a  simu¬ 
lation. 

The  numerical  results  for  the  case  of  M=2  with  X^  ■  X^  and  p^  »  p^  are  shown  in 
Fig.  4  and  with  X^  **  Pj  *  P2  in  Fig.  5.  Numerically  optimizing  over  the  values 
of  the  p^'s  we  find  that  in  the  non-symmetric  loading  case  the  delay  performance  is 
consistently  superior  over  all  throughput  values  to  that  of  the  best  symmetric  case, 
see  Fig.  6.  In  fact  the  optimum  non-symmetric  case  requires  that  one  of^the  two  p^s 
be  equal  to  one.  The  other  p^-value  can  be  computed  and  varies  with  X.  The  implica¬ 
tion  is  that  by  favoring  one  of  the  two  users  the  average  delay  is  improved.  Thus 
the  gains  achieved  by  one  user  outweigh  the  losses  incurred  by  the  other.  This 
observation  may  be  of  value  in  the  design  of  access  protocols  if  it  is  explained  and 

confirmed  by  rigorous  theoretical  justification. 
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To  evaluate  this  approximate  model  a  Monte-Carlo  simulation  was  run  for  a  2-node 
and  a  3-node  system  with  various  values  of  retransmission  probabilities  and  arrival 
rates.  The  comparison  to  the  computational  results  is  shown  in  Figures  7  and  8 

A 

respectively,  for  non-symmetric  loading  cases.  The  average  delay  for  each  terminal 
1 8  plotted  separately.  The  arrival  rates  are  equal  for  all  terminals. 

In  Fig.  9  we  compare  the  approximate  results  from  our  model  to  the  analytical 
results  of  [a].  The  case  of  two  users  with  symmetric  loading  *s  t*ie  only 

one  that  has  been  solved  exactly  and  thus  such  a  comparison  is  necessary.  We  find 
excellent  agreement,  with  difference  between  the  two  sets  of  results  bounded  by  1.51. 
Equally  good  agreement  is  found  between  optimal  delays  for  the  approximateand  analy¬ 
tical  cases  for  the  symmetric  two-user  system  as  shown  in  Fig.  10. 

7 .  Conclusion 

In  this  paper  a  variation  of  an  approximate  model  for  the  delay  analysis  of  a 
slotted-ALOHA  system  with  a  finite  number  of  buffered  terminals  was  proposed  and  ana¬ 
lyzed.  Limited  performance  evaluation  was  carried  out  for  small  numbers  of  terminals 
in  terms  of  comparison  of  the  results  to  the  analytically  obtained  ones  (whenever 
exact  analysis  is  available)  and  to  those  obtained  by  simulation.  Agreement  proved 
to  be  very  good.  Not  having  evaluated  the  model  for  greater  values  of  numbers  of 
terminals  we  can  only  be  moderately  confident  in  its  predictive  capability.  Our  con¬ 
fidence,  however,  is  backed  by  the  arguments  provided  about  the  structural  nature  of 
the  model  that  served  also  as  our  primary  source  of  motivation  in  proposing  it.  We 
believe  it  may  prove  to  be  a  useful  and  usable  tool  in  evaluating  performance  of 
several  protocols  (such  as  CSMA)  that  rely  on  contention  among  a  finite  number  of 
terminals  that  maintain  packet  queues. 

*  the  little  "square"  and  "triangle"  marks  on  the  curves  do  not  represent  the  points 
on  the  basis  of  which  the  curves  were  drawn  but  rather  serve  to  identify  the  curves. 
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APPENDIX  1 


TRANSITION  PROBABILITIES  OF  SYSTEM  STATUS  CHAIN 
In  this  Appendix  we  calculatethe  transition  probabilities  of  the  system  status 
chain  as  discussed  in  section  3.1.  We  assume  that  the  quantities  s^  and  r^  are  equal 
to  one;  that  is,  a  packet  gets  immediately  transmitted  when  the  terminal  is  idle  or 
active.  The  same  calculations  can  be  easily  carried  out  when  this  is  not  the  case. 
Also  note  that  p  and  A  stand  for  1-p  and  1-A  respectively. 

Since  the  aim  is  to  develop  the  steady  state  distribution  for  the  system  status 
chain,  we  need  to  calculate  the  transition  probabilities.  These  calculations  reduce 
to  considering  a  small  number  of  cases,  based  on  the  change  in  the  number  of  active 
nodes,  thus  we  classify  the  transitions  into  four  types  as  follows: 

1)  The  number  of  active  nodes  transits  from  0  to  1 

An  idle  node  cannot  transit  into  active  status,  so  the  only  possibility  for  such 
a  transition  is  one  involving  a  blocked  node,  say  node  j,  which  becomes  active,  while 
the  other  blocked  nodes  do  not  attempt  to  retransmit  and  all  idle  nodes  receive  no 
packets.  The  probability  of  such  a  transition  is  expressed  by: 

M-n-1  _  n+1  _  p, 

P(AA  -  1,AB  -  -  1)  -  n  A  n  p.  =J-  [(1  -  P.(0|2))  +  A  P  (0|2)  ] 

k-1  K  i-1  1  Pj  3  3  J 

M-n-1  _  n+1  _  p,  _ 

-  "All  p  =■*-  [l-P  (0 12))  A  ]  (1.1) 

k-1  i-1  Pj  3  3 

where  AA  and  AB  represent  the  changes  in  the  numbers  of  active  and  blocked  nodes 

and  the  first  term  in  brackets  is  the  probability  that  the  buffer  content  of  the 

blocked  node  j  is  not  0,  while  the  second  term  is  the  probability  that  the  buffer 

content  of  node  j  is  0  and  a  packet  arrives.  All  other  quantities  Involved  have  been 

defined  in  section  3.  Note  that  for  ease  of  notation  the  subscripts  i  and  k  are 

shown  to  vary  from  1  to  n+1  and  from  1  to  M-n-1  while  it  is  meant  that  they  vary  over 

the  sets  of  blocked  and  idle  nodes  respectively.  The  same  notation  is  followed  in 
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the  sequel 


2)  The  nunber  of  active  nodes  remains  1 

In  this  case,  the  status  of  all  Idle  nodes  and  blocked  nodes  must  be  kept 
unchanged.  The  probability  of  such  a  transition  is: 

M-n-1  _  n  _ 

p(aa-o)-  n  x  n  p  [(l-p  (i|D)  +  x,p.(i|D]  •  p, 
k-1  i-1  2  2  2  2 

M-n-1  _  n  _ 

-  n  X  n  p± [1-P,(1 |i) X  ]p  (1.2) 

k=l  i=l  J  J  J 

where  the  first  term  in  the  brackets  is  the  probability  that  the  buffer  content  of 
the  active  node,  node  j,  is  greater  than  one  and  the  second  term  is  the  probability 
that  the  buffer  content  is  one  while  one  packet  arrives. 

3)  The  nunber  of  active  nodes  transits  from  1  to  0 

We  subdivide  this  sort  of  transition  into  three  sub-types.  The  first  sub-type 
is  that  the  active  node  turns  to  idle.  The  probability  of  such  a  transition  is  given 
by: 

M-n- 1_  n 

P(AA  -  -  1,  AI  -  1)  -  n  ^  H  Pl[l-P  (1  |l)  Jx  ,p  (1.3) 

k=l  i=l  J  J  J 

The  second  sub-type  is  that  the  active  node  packet  collides  with  packets  from 
blocked  node(s)  and  becomes  blocked.  The  probability  of  the  transition  is: 

M-n  _  n-1 

P(AA  -  -  1,  AB  -  1)  -  n  X  (1  -  n  p  )  *  p  (1.4) 

k=l  i=l  2 

The  third  sub-type  is  that  the  active  node  collides  with  q-1  idle  nodes,  which 
results  in  the  active  node  and  the  q-1  idle  nodes  turning  to  blocked.  The  probabi¬ 
lity  is: 

q- 1  M-n  _ 

P(AA  -  -  1,  AI  -  -  (q-l),AB  -  q)  -  n  II  •  p  (1.5) 

kj-1  1  k=q  2 
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4)  The  number  of  active  nodes  remains  0 


In  this  case  there  are  four  sub-types.  Hie  first  one  is  that  all  nodes  maintain 
their  status  without  change.  Hie  transition  proability  is: 

n  _  n  p.  M-n  _  n  _  M-n  _  M-n  X, 

p(ai  -  o,  ab  -  o)  -  [1  -  n  p  z  (-J-)  ]  nx+np  UXE  (•*£) 

i=  1  1  j-1  Pj  k=l  *  i-1  1  k-1  k-1  \ 

M-n  _  n  _  n  p  M-n  X 

+  n  X  [l  -  n  p  (  Z  (=1)  -  Z  Or^))]  (1.6) 

k=l  K  i=l  j«=l  pj  k=l  \ 

where  the  first  term  is  the  probability  that  no  packet  arrives  at  the  idle  nodes 
while  no  blocked  node,  or  at  least  two  blocked  nodes,  transmit,  and  the  second  term 
is  the  probability  that  no  blocked  node  attempts  to  transmit  while  only  one  of  the 
idle  nodes  receives  an  arriving  packet. 

The  second  sub-type  is  that  one  of  the  blocked  nodes,  say  node  j,  becomes  an 
idle  node.  The  transition  probability  is: 


M-n-1  __  n+1  _  p, 

P(AI  -  1,  AB  -  -  1)  -  n  X  n  p.  [X,  P,(0|2)  ]  (1.7) 

k=l  K  i=l  1  3  Pj  3 


where  the  quantity  in  the  brackets  is  the  probability  that  the  blocked  node  attempts 
to  transmit  while  its  buffer  content  is  0  and  no  packet  arrives  at  the  node. 

The  third  one  is  that  only  one  of  the  idle  nodes,  say  node  j,  becomes  a  blocked 
node.  The  probability  is: 


P(  AI 


AB  =  1) 


M-n-1 

n 

k=l 


X.  n-1  _ 

\^(1  -i-ni  ''i) 


(1.8) 


Hie  fourth  one  is  that  q  (q>2)  of  the  idle  nodes  become  blocked.  Hie  transition 
probability  is: 


q  X.  M-n  _ 

P(AI  -  -  1,  AB  -  1)  -  n  (^2-)  .n  X  M-n>l  (1.9) 

j-1  j 
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Based  on  these  transition  probabilities  we  can  express  explicitly  the  state 
equations  for  the  system  status  state.  Note  that  they  will  depend  on  known  quan¬ 
tities  (X^'s  and  p^'s)  and  the  unknown  quantities  P^(0|2)  and  P^(l|l).  In  the  next 
appendix  these  quantities  will  be  expressed  in  terms  of  the  state  probabilities  of 
the  node  chain,  ir(T^,N^). 
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APPENDIX  2 


PROBABILITY  DISTRIBUTION  OF  QUEUE  LENGTH 

In  this  Appendix,  we  prove  equations  (5)  -  (10)  according  to  the  state  tran¬ 
sition  diagram  of  Figure  11.  Here,  too,  we  assume  *  r^  m  1,  and  we  denote  *  p^. 


The  steady  state  equations  for  node  i  are  as  follows: 

it(0,N)  “  \PB(i)»(0,N)  +"XiPA(i)ir(l,Nfl)  +  1^(1)  ir(l  ,N) 

+  X^Pg(i)  ir(0,N-l) ,  N>0  (2.1) 

it(l,N)  =  A^UM'l.N)  +\PB(i)Tr(0,N)  +  A^D^O.N-l) 

+  N>0  (2.2) 

tr(0,0)  -TJB(i)«(0,0)  +‘AiPA(i)r(l,l)  +  A^i)  *(1 ,0)  (2.3) 

x(l,0)  -  (1  -  A^?1(i)  )x(l  ,0)  +"AiPB(i)ir(0,0)  +  7^(1)  ir(l  ,1)  (2.4) 


From  the  equations  of  ir(l,0)  and  ir(0,0),  we  obtain  the  following  equations 
after  multiplication  of  both  sides  by  z^  and  addition: 

[(PB(i)  +  A±  -  Pfi(i)Ai)  -  A±(l  -  Pfi(i)  )z  ]Gq(z)  - 

[ AiPB(i)  +“AiPA(i)z"1]Gj(z)  +  [A±(PA(i)  -  PI(i))  -  ]<1 ,0)  (2.5) 

From  the  equations  for  ir(l,0),  ir(l,n),  we  obtain 

[l  -  A^Fg(i)  -l^Dz'^G^z)  -  r^PgU)  +  AiPg(i)z  ]Gq(z) 

+  [l  -  AjjU)  -  \±  +  A^i)  -  A^Dz^Jnd.O)  (2.6) 

Eq.  (2.5)  yields  for  z  «  1 

PB(i)G^(l)  -  Pa(1)gJ(1)  +  (AjPjU)  -  PA(i))x(l,0)  (2.7) 


20 


After  differentiating  (2.5)  we  obtain  at  z  ■  1: 


PB(i)GQ(l)  +  VB(i)Go(l)  "  Vi)Gl(1)  "  XiPA(i)Gl(1)  "  "  XiPA(1)l,(1’0)  (2,8) 

By  differentiating  (2.6),  evaluating  at  z  ■  1,  and  adding  to  (2.8),  we  obtain 


Gl(1)  “  Xi  +  X1ir(1»0) 


(2.9) 


Then  from  (2.7),  (2.9),  and  the  fact  that  Gq(1)  +  G*(l)  *  1  we  obtain  (6),  (8), 
(9).  From  the  equations  for  ir(0,0)  and  tt(1,0),  we  obtain  (5),  (7).  We  use  the 
generating  property  of  G  to  calculate  the  average  length  L^»  namely 


Li  "  Go(1)  +  Gl(1) 


(2.10) 


In  order  to  calculate  these  derivatives,  we  use  (2.8),  which,  after  substitution 
from  (5)  to  (9)  for  values  of  ir(0,0),  it(1,0),  ir(l,l),  G^(l)  and  gJ(1),  becomes 


M  X.X2PT(i)(P_(i)  -  P.(i)) 

P.(i)Gn(l)  -  PA(i)G,(l)  =  - - - - - 

X  P  (i)  -  X  (P  (i)  -  P  (i)) 


(2.11) 


Also,  we  add  (2.5)  and  (2.6),  differentiate  twice,  and  evaluate  at  z 
we  obtain 

X. 


1.  Thus 


gJcu 


•i  -i  X  X  P  (i) 

-jr-Gjd)  +- - - 

Ai  1  X  P  ( i)  -  X  (P  (i)  -  Pt(i)) 


i  B' 


i'  I' 


(2.12) 


From  (2.11),  (2.12)  we  have 


xfri^rw1) +  xipa(i)^ 


Gjd) 


gJ(1) 


(X^i)  -  X±PA(i))f^PB(i)  -  X1(PI(i)  -  PAd))] 

x?ipi(i)  rxiPB(i)  +  xipA(i)^ 

(X^jjd)  -  X^PA(i))  [^Pg(i)  -  X1(PI(i)  -  PA(i))] 


(2.13) 


(2.14) 
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Thus 


L±  *  Gj(l)  +  G?(l) 


XWlW 


1  <VB(1)  "  Va(1^Vb(±)  "  Xi(PI(±)  "  PA(i))^ 


(2.15) 


which  is  Eq.  (10). 

Finally  from  the  definitions  In  section  3.1  it  follows  that 


PjU  |D 

P±(0 |2) 


*(1,1) 


Gj(l)  -  *(1,0) 


Tr(O.O) 

gJ(1) 


By  using  (5)  -  (7) ,  we  prove  also  that 

Pi(l  |l)  -  P1(0|2)  (2.16) 

Thus  all  equations  now  are  explicitly  obtained  and  permit  the  application  of  the 
iterative  algorithm  shown  in  Fig.  3. 
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Figure  Is  Model  of  a  Slotted  ALOHA  System 


Figure  2  ;  System  Status  Markov  Chain  for  M*2;  each  pair 
denotes  the  status  of  the  two  terminals;  0,1,2 
stand  for  idle,  active,  and  blocked  respectively. 


Figure  3  :  Iterative  Algorithm  for  the  solution  of 

the  coupled  sets  of  equations  for  the  state 
probabilities  of  the  system  chain  and  the 
queue  length  chains. 
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Figure  4  ;  System  Delay  vs .  Throughput  for  Symmetric  Case,  M*2 
curves  are  parametrized  by  p(Pj~p2~^* 
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Figure  5  •  Average  System  Delay  vs.  Throughput  for  Non-symmetric  Case,  M=2 
Curves  are  parametized  by  the  pairs  of  values  of  p^  and  p2 
respectively. 
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Figure  6  :  Optimal  Delay  vs.  Throughpu 
in  Symmetric  and  Non-Symmet 
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Figure  7  ;  Comparison  Between  Numerical  Results  and  Simulation 


(M  «  2) 
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Figure  9  Delay  vs.  Throughput  for  Symmetric  Case  of 
Scheme  2  (M  «  2) 
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Figure  10:  Comparison  Between  Optimal  Delay  vs.  Throughout 
Between  Numerical  and  Analytical  Results 


Figure  il  :  Transition  of  Node  Status  and  Queue  Length 


